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1. INTRODUCTION 

The depletion of fossil fuels and the increasing energy demand over the world, about 55% by 
2035 [1], are brought more attention to green energy. The production of electricity from conventional sources 
affects environment balance and causes pollution. This pollution affects lifes and animals. The power produced 
from fossil fuels gives off harmful gazes such as oxides of carbon. These gazes contribute to global warming. 
Limited reserves of fuels and their unstable costs are the most important reasons for renewable energy. 
The increasing concern about environmental pollution and the impact of traditional sources has emphasizes all 
countries for reducing their emission. Renewable energy sources seem to be the best solution for a sustainable 
electrification. We can have the electricity directly from sunlight via PV panel. Solar energy is one of the most 
promising renewable energy technologies: it is clean and abundance. However, the intermittent nature of these 
types of energy makes the energy produced from one renewable source unreliable [2, 3]. Coupling PV panel 
with another source of energy such as wind turbine generator can reduce significantly the intermittence issue. 

The electric Hybrid renewable energy system is less costly and more reliable than system with one 
source [4, 5]. Designing a hybrid renewable energy system is a difficult task, the sizing or the number of 
element and the used control strategy is very essential, the performance of the PV/Wind/Fuel cell system can 
be significantly influenced by the proposed control strategy. Wind turbine generator and PV panel might 
operate in the way that they can give the maximum power to achieve a high efficiency value [6]. In this study, 
the photovoltaic (PV) panel and the wind turbine generator are used as the main source for the load demand. 
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the fuel cell (FC) electrolyzer is required as a storage capacity. The photovoltaic panel (PV) technology is seen 
as the easiest renewable source due to several reasons: The study of a hybrid renewable energy system (HRES) 
coupled with hydrogen energy is became an area of interest [7]. 

The traditional storage with chemical batteries provides high discharging efficiency and high energy 
storage capacity [8]. This system provides an improvement to overcome some drawbacks of RES especially 
the intermittence problem: the energy storage systems are able to feed load demand or consume the produced 
energy from renewable sources. Diesel generators or fuel cells seems necessary in hybrid renewable energy 
systems (HRES) by supplying the load demand when the storage facilities are enable or empty. The paper 
presents IBSA (improved backtracking search algorithm) to optimize the annual total cost (ATC) for 
photovoltaic generator, wind turbine generator and fuel cell system. 


2. HYBRID RENEWABLE ENERGY SYSTEM PRESENTATION 

The system used in this paper is a multisource system composed of hybrid renewable energy source 
system based on photovoltaic panel, wind turbine generator and fuel cells electrolyzer [9]. The storage facilities 
are used to smooth out the renewable energy fluctuation. Besides, storing electricity in large scale during 
off-peak hours reduce significantly the dependence of fossil energy during peak demand. 


2.1. Photovoltaic (PV) panel model 

To determine the power output for PV panel, two models are possible: the probabilistic one and 
the deterministic one. Both are based on climatic data. In this study we use the deterministic model to define 
the PV power. The electric power is determined directly from the nominal power of the PV cell. This power is 
injected directly in the DC bus: 


Faire = Eair ` COS 0 (1) 


where Edirc is the direct sunlight received by the collector, Ois the inclination angle. 





1+cos 0 
Eqife = Eaif ` (=) (2) 
Edifc is the diffuse sunlight received by the collector. 
Ppvc 
Poy = (Eire + Eaifc) i F (3) 
Ppvc is the nominal power of PV panel and En is the insolation in standard conditions. 
cos 0 = cos $ : cos(os, gc) + sin È» sin $ : cos È (4) 


B is the altitude of the sun, sis the azimuth, cis the collector azimuth and > is the inclination angle of 
the solar collector. 


sin p = cos L » cos ô : cos H + sin L » sin ô (5) 
L is the latitude of the location, ò is the solar declination and H is the hour angle. 


cos 6:‘sinH 


(6) 


sin pS = oe 
2.2. Wind turbine model 

Wind turbines use the kinetic energy of wind speed. In this paper the power delivered by wind turbine 
generator is modeled according to wind speed, as we can see in Figure 1. The power output of 
the wind turbine generator is based on wind speed. It is expressed by (7) as we can see in (7): 


0, V < v4 
1 2 3 
—--o:17:R*“:-V°?, vaxsv<v 
W = R d n (7) 
Po a a SV, 
0 v > Ve 


where Vd is the boot speed, Vn is the nominal speed and Vc is the shut down speed. 
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Figure 1. Wind turbine output power 


2.3. Fuel cell model 

There are different types of Fuel cell (FC); the most used for distributed generation is the proton 
exchange membrane (PEM) [10]. In this work, we consider PEM FC; it is studied by using the model of 
the FC stack. The output voltage of the FC Vfc can be calculated from the voltage developed inside the FC. It 
is expressed by the (8): 


Vee = Ege m Vat Vohm (8) 


where Eoc is the voltage when the circuit is open, Vact is the intern voltage of the FC and Vohm is the ohmic 
voltage. 
The number of hydrogen storage tanks is determined based on the Management strategy 
showed below: 
— if the power generated from renewable energy sources (PV/Wind) is greater than the load demand, 
the electrolyser will be used to produce hydrogen. This hydrogen stored in the tanks can be calculated 
as follows: 


E stor (t) = E stor (t 7 1) + (Eren (t) — Eoad (E)/EF finv) Ef fetect (9) 


where Estor is the hydrogen stored in the tanks, Eren is the energy produced from renewable sources, Eload 
is the energy load and Effiny is the efficiency of the inverter. 

— when the energy demand of the load is greater than the energy produced by renewable energy sources, 
the FC will be used to feed up the load [11]. The amount of hydrogen stored in the tanks is calculated 
as (10). 


Estor (t) = Estor (t — 1) — (Etoaa(t) — Einv(t))/Ef fruet ceu (10) 
The cost of supplying the energy demand by FC electrolyzer is determined by the (11): 


CFC 


eee + Cogm Fc (11) 


Csrc = 


where Crcis the FC purchase cost, LifeFC is the FC lifetime and Cogm rc are the FC operation and maintenance 
costs. 


3. BACK-TRACKING SEARCH OPTIMIZATION ALGORITHM (BSA) 

BSA is a methaheuristic algorithm, based on natural or biological evolutions techniques such as 
mutation. It is proposed to solve constrained optimization problems and overcome some drawbacks of 
the previous evolutionary algorithms: e.g., high sensitivity to the control parameter and time-consuming 
computation [12]. BSA’s structure is simple; it has a powerful global exploration and local exploitation due to 


TELKOMNIKA Telecommun Comput El Control, Vol. 18, No. 1, February 2020: 456 - 464 


TELKOMNIKA Telecommun Comput El Control O 459 


controlling the search direction by scale factor parameter [13, 14]. Figures 2 (a) and (b) shows 
power as a function of current and the voltage as a function of current for the FC used in this study. Figure 3 
shows the general structure of the BSA algorithm. BSA has five main steps: initialization, selection-I, mutation, 
crossover and selection-II [15]. Its general structure is described in Figure 3. Figure 4 shows the initialization 
process of BSA. 
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Figure 2. (a) Power-current (b) Voltage-current curves of the FC 


Initialization Initialization 
Repeat if numel(low)==1, low=low*ones (1,dim); 
Selection-I 
Generation of trial population up=up*ones (1, dim) ; 


Mutation 


Crossover end 
End pop=GeneratePopulation (popsize,dim, low, up) ; 
Selection-II fitnesspop=feval (fnc,pop,mydata) ; 
Until stopping criteria are met historical pop=GeneratePopulation (popsize, dim, low, up) ; 








Figure 3. General structure of BSA Figure 4. BSA’s initialization 


— Initialization 
In this process BSA generates the initial population by a uniform distribution. 
— Selection-I 


BSA defines the Old population Old P, which is used to calculate the search direction. The historical 
population is initialized by (12). 


In each iteration, BSA determine Old P by comparing two numbers generated randomly a and b. This 
comparison is shown in (13): 


if a< b then oldP = P\a, b/U(0,1)| (13) 


after determining the historical population Old P, BSA changes the order of individuals by using 
the permute function. 
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oldP = Permutting(oldP) (14) 


— Mutation 
In this step BSA generates the initial form of trial population. 


Mutant = P + F * (old P — P) (15) 
Where F is a function. Its value controls the amplitude of the search direction. It is expressed by (16): 

F = 5 *randn, where randn = N(0,1) (16) 
— Crossover 


M apirana(p) =1 

M ap iuc1:[-mixrate+rand«D]) = 1; u = permuting (1,2,3,..., D) 

Crossover = P + (map. F).* (old P — P) 

Tij = rand * (up ie low; ) + low, if T; ; beyond the boundary. (17) 


— Selection 
In this step, the individual with best fitness value is used to replace the previous individual. 


4. IMPROVED BACKTRACKING SEARCH OPTIMIZATION ALGORITHM (IBSA) 

BSA has a random strategy in defining mutation; the later introduces changes in individual position. 
The mix rate parameter (M) used in crossover operator controls the number of individuals that mutate in trial 
population. In this work the improved backtracking search optimization algorithm (IBSA) presented by [16] is 
used to study the impact of using a FC electrolyzer in a hybrid renewable energy system especially in terms of 
energy cost. 

The IBSA is proposed to overcome some drawbacks of the traditional BSA such as the convergence 
issue [17-19]. However, BSA shows a strong robustness in finding the best cost value; its way in storing 
population from the previous generation may make it converge very slowly [20, 21]. The system used in this 
study is presented in Figure 5. It’s composed of PV generator, wind turbine generator and Fuel cell system. 
The idea behind the IBSA is defining a new mutant based on the scale factor value. To see 
the effectiveness of the IBSA four benchmark functions are used as shown in Table |. The simulation results 
are shown in Figure 6 and the statistics values are summarized in Table 2. 





Figure 5. PV wind turbine generator fuel cell hybrid renewable system 
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Table 1. Low-up and dimension for the benchmark functions 





ID Name Low Up D 
Fl Ackley -32 20.39 20,18 
F2 Criewank -600 0.2423 0.3449 
F3 Schaffer -100 20.24 19.63 
F4 Shekel -32 1.824 1.643 
OF 
-0.002 
-0.004 | 4 
: -0.006 # 
= = 0.008 
8 8 
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(c) (d) 
Figure 6. Optimal cost value for (a) F1, (b) F2, (c) F3 and (d) F4 using 
BSA (green curve) and IBSA (red curve) 


4.1. Objective function 
The ATC of the system can be defined as follows: 


y1 = C8. Pay + C4.W + CE. Pre (18) 


where Cp is the total cost of PV generator, Cw is the total cost of wind turbine generator and Cfc is the total 
cost of the fuel cell. The cost of supplying energy with PV generator and wind turbine generator is determined 
in [22]. The cost of power produced from FC is determined in (19). 


CFC 
Liferc 





Csrc = + Cogm_Fc (19) 


4.2. Electrolyzer model 
The hydrogen production nh2 rate is determined based on Faraday’s law [23], it is proportionnel to 
the electrical current inside the circuit [24]. The nh2 (mol/s) is determined in (20). 


Nya = “Ee (20) 


where np is the Faraday efficiency, nc is the number of cells in series ieis the electrolyzer current and F is 
the Faraday constant. 
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5. SIMULATION AND RESULTS 
The studied system is composed of PV generator, wind turbine generator and FC electrolyzer. 
The load demand is an AC load; and the inverter is used with an efficiency of 0.8. The FC is used to feed up 
the load in worse case for RE (renewable energy). The system flowchart is presented in Figure 5. The used 
system is composed of [25]: 
— 6000 Wp PV generator, with a total acquisition cost of 40000 £. 
— 5000 W wind turbine generator, with a total acquisition cost of 10000£ and an annual operation and 
maintenance cost of 300£. 
— IKW FC, acquisition cost of 4000£, O&M cost of 0.2£/h, 30000 h expected lifetime, 
Nfc = 16.66 KWh/KgH2 for the nominal power and Pmin = 60 W. 
— IKW electrolyzer, acquisition cost of 3200£, Neyz = 0.021 KgH2/KWh. 
— 0.1Kg H2 tank, acquisition cost of 150 £, O&M cost 10£/year, 25 years expected lifetime. 
The main objective of this work is to size optimally a hybrid PV/Wind/FC system using IBSA and 
compare the result with the BSA. The proposed method aims to satisfy many requirements: 
— Optimize the annual total cost (ATC) of the system. 
— Demonstrate the effectiveness of the IBSA compared to the traditional BSA for system using Fuel cell. 
— Study the influence of the electrolyzer current on the Faraday efficiency. 
Table 1 and Table 2 show respectively the statistics values and the low-up and dimension for 
the benchmark functions. 


Table 2. Statistics values for benchmark functions 


Problem Statistics BSA IBSA 
Fl Mean 20.39 20.18 
Std 0.2423 0.3449 

Best 20.24 19.63 

F2 Mean 1.824 1.643 
Std 0.1014 0.093 

Best 1.753 1.625 

F3 Mean 0.16 0.067 
Std 0.0573 0.1168 

Best 0.1241 0.0048 

F4 Mean -0.010 -0.009 
Std 0.0041 0.002 

Best -0.015 -0.011 


5.1. Annual total cost optimization of hybrid system 

Figure 5 shows the configuration of the Hybrid renewable PV/Wind/FC system under study. 
The proposed structure meets the power demand for isolated load located in Rabat, Morocco. Figure 6 plots 
the optimal value for the four benchmark functions for BSA and IBSA. In this work the PV generator and wind 
turbine generator are the primary energy source, while fuel cell electrolyzer is used as a backup system. 
The optimization method based on IBSA (improved backtracking search optimization algorithm) is used to 
size and design optimally the hybrid renewable energy system. The program is developed on MATLAB 
software. IBSA optimizes the numbers of PV panels (Npv), the numbers of wind turbines (Nw) and 
the numbers of FC electrolyzer (Nfc). These numbers are included in the power value of each component as 
shown in (18). 

The IBSA must meet the load demand with a minimal cost value. The optimal annual total cost (ATC) 
is shown in Figure 7. The curve with green color is for BSA and the one with red color is for IBSA. To make 
a good comparison, we have maintained the same parameters for BSA and IBSA, 1.e., population size and 
the number of generation. The simulation process is initialized by a random value of Npv, Nw and Nfc. It is 
demonstrated that the best ATC obtained by IBSA is better than the value obtained by BSA. 

This work shows also the impact of the electrolyzer current on the Faraday efficiency Figure 8. 
The Faraday efficiency increases with the augmentation of the current that is means a high value for 
electrolyzer current gives a high level of efficiency. On the other hand, the renewable system has designed to 
feed up the load demand. This can be achieved through the electrolyzer current whish is directly linked to 
the produced hydrogen. 

The level of hydrogen stored in the tank influences the optimization process (the value of ATC) and 
then the number of PV panels, the number of wind turbine and the number of FC. The number of PV panels 
Npv, the number of wind turbine generator Nw and the number of FC Nfc optimized by IBSA is shown in 
Figure 9. The optimal sizing and costs of the Hybrid PV/Wind/FC system is as follow: 
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— 6KWof PV 

— 5 KW of Wind 

— 1 KWof FC 

with an Annual Total Costs (ATC) of 14562.17 £. 


ATC(E) 





0 20 40 60 80 100 120 140 160 180 200 
Iteration 


Figure 7. Optimal ATC with BSA (curve green) and IBSA (curve red) 
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Figure 8. Faraday efficiency-current curve Figure 9. PV number/wind turbine/FC 


6. CONCLUSION 

This work has used a novel optimization method based on Improved backtracking search optimization 
algorithm (BSA). This is the first time when this new algorithm is used to perform a technical optimization for 
system compromising PV panels, wind turbine generators and fuel cell (FC). Firstly, a comparative study is 
done based on four benchmark functions: ackley, criewank, schaffer and schekel. The result shows 
the effectiveness of the IBSA and Table 1 summarizes the statistical value for the two compared algorithm 
BSA and IBSA. IBSA outperforms BSA in terms of convergence speed and best fitness value. Secondly, 
the system under study is modeled through several power equations then the objective function (OF) is defined. 
We have optimized the OF using BSA and IBSA, and it can be easily seen that the IBSA gives the best ATC 
value. Finally, we have concluded that the Faraday efficiency influences the electrolyzer current and so 
the hydrogen stored in the tank. 
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